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Abstract 



We study the properties of static and rotating neutron stars adopting 
non-relativistic equations of state (EOS) for asymmetric nuclear matter 
based on the Brueckner-Hartree-Fock (BHF) scheme. The BHF calcula- 
l/"") | tion, with the continuous choice for the single particle potential, appears 

to be very close to the full EOS, which includes the three-hole line con- 
■ tribution calculated by solving the Bethe-Fadeev equations within the 

gap choice for the single particle potential. Three-body forces are in- 
cluded in order to reproduce the correct saturation point for nuclear 
matter. A comparison with fully relativistic many-body calculations of 
nuclear matter EOS is made. 
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3 ■ 1 Introduction 

H ' 

The properties of neutron stars such as masses and radii depend on the equa- 
tion of state (EOS) at densities up to an order of magnitude higher than those 
observed in ordinary nuclei. Therefore the knowledge of the EOS in the su- 
perdense regime is fundamental for the study of astrophysical compact objects 
[|]]. For this purpose we derive an equation of state for asymmetric nuclear 
matter, using a non-relativistic many-body theory within the framework of 
the Brueckner-Hartree-Fock (BHF) scheme @, §. In this approach, the basic 
input is the two-body nucleon-nucleon (NN) interaction. The BHF approxi- 
mation, with the continuous choice for the single particle potential, reproduces 
closely the many-body calculations up to three hole- line level || . However, as 
it is well known, the empirical saturation point is not reproduced. Therefore 
we have included a contribution coming from three-body forces to reproduce 
the correct saturation point ||. Those EOS's are the fundamental input for 
constructing models of static and rotating neutron stars in the framework of 



Einstein's theory of general relativity by applying a refined version of Hartle's 
stellar structure equations || 0. We calculate properties of neutron stars like 
gravitational mass, equatorial and polar radius, for sequences of star models 
either static or rotating at their respective general relativistic Kepler frequen- 
cies. We compare with predictions from fully relativistic microscopic EOS. 



2 Equation of state 

Microscopic calculations of nuclear matter EOS have been performed in the 
framework of the Brueckner-Hartree-Fock (BHF) scheme || . The energy per 
particle E/A within the BHF scheme is given in terms of the so-called reaction 
matrix G. The latter is obtained by solving the Brueckner-Bethe-Goldstone 
(BBG) equation 

G(u) = V + V — Q—G(u), (1) 

U) — Hq 

where u is the unperturbed energy of the interacting nucleons, V is the free 
nucleon-nucleon (NN) interaction, H is the unperturbed energy of the inter- 
mediate scattering states and Q is the Pauli operator which prevents scattering 
into occupied states. With the G-matrix we can calculate the total energy per 
nucleon 

A 5 2m v ; ' v ; 

being U(n) the contribution of the potential energy to the total energy per 
particle 

U ^ n ) = VA £ (kk'\G(uj = e(k) + e(k'))\kk') a (3) 
/yi k,k'<k F 

where the subscript a indicates antisymmetrization of the matrix element. 
The single-particle energies are denoted by e. In this scheme, the only input 
quantity we need is the bare NN interaction V in the Bethe-Goldstone equation 

(I)- 

The Brueckner-Hartree-Fock (BHF) approximation for the EOS in sym- 
metric nuclear matter, within the continuous choice H, reproduces closely 
results which include up to three hole-line diagram contributions to the BBG 
expansion of the energy, calculated within the so called gap choice for the 
single particle potential 0]. In Fig.l we show the energy per nucleon calcu- 
lated within this scheme in the case of symmetric and neutron matter using 



the Argonne v±4 model || for the two-body nuclear force. The open squares 
represent the solution of the Bethe-Fadeev equations for the three hole-line 
in the gap choice. These results extend to higher densities the calculations 
published in ref.0] for symmetric nuclear matter. Recently it has been shown 
that the results up to three hole-lines are independent from the choice of the 
single particle potential ||, which gives evidence of convergence of the BBG 
expansion. 

For neutron matter the results are preliminary. In any case we found a 
negligible difference between the BHF results in the continuous and in the gap 
choice. Correspondingly the three hole-line contribution in neutron matter 
turns out to be much smaller than in symmetric nuclear matter. All these 
results together give support to the use of the BHF approximation in the 
continuous choice in the study of neutron stars. 

We notice that the BHF fails to reproduce the empirical saturation point 
of nuclear matter. This well known deficiency, which does not depend on 
the choice of the two-body force, is commonly corrected introducing three- 



body forces (TBF). We adopted the Urbana three-nucleon model [10|, which 



consists of an attractive term due to two-pion exchange with excitation of 
an intermediate A-resonance, and a repulsive phenomenological central term. 
Several details are given in ref.f|. 

The corresponding EOS obtained using the t>i 4 potential is depicted in 
Fig. 2a) for symmetric and neutron matter (solid line). This EOS saturates at 
n Q = 0.178 fm~ 3 , E /A = —16.46 MeV and it is characterized by an incom- 
pressibility K^, = 253 MeV. In Fig. 2a) we plot also the EOS from a recent 
Dirac-Brueckner calculation (DBHF) |TT| with the Bonn-A two-body force 
(dashed line). In the low density region the BHF equation of state with TBF 
and DBHF equation of state are very similar, whereas at higher density the 
DBHF is stiffer. The discrepancy between the non-relativistic and relativistic 
calculation of the EOS can be easily understood by noticing that the DBHF 
treatment is equivalent to introduce in the non-relativistic BHF the three- 
body force corresponding to the excitation of a nucleon-antinucleon pair, the 
so-called Z-diagram which is repulsive at all densities. In BHF treatment, on 
the contrary, both attractive and repulsive three-body forces are introduced, 
and therefore a softer EOS is expected. 

The properties of neutron stars (NS) depend on the knowledge of the EOS 
over a wide range of densities, i.e. from the density of iron at the star's 
surface up to several times the density of normal nuclear matter encountered 
in the core |JJ. It is commonly accepted that the interior part of a neutron 
star is made mainly by nuclear matter (eventually superfluid) with a certain 



lepton fraction, although the high-density core might suffer a transition to 
other hadronic components. Here we assume that a neutron star is composed 
only by nucleons and leptons, i. e. an uncharged mixture of neutrons, protons, 
electrons and muons in /^-equilibrium. The presence of leptons softens the EOS 
with respect to the pure neutron matter case. In Fig. 2b) the EOS for /3-stable 
matter is shown for the cases previously discussed. It has to be stressed that 
in /3-stable matter the density dependence of the nuclear symmetry energy 
affects the proton concentration ||13||. As it has been recently pointed out by 



Lattimer et al. ||14|| , the value of the proton fraction in the core of NS is crucial 
for the onset of direct Urea processes, whose occurrence enhances neutron star 
cooling rates. 

From the energy per baryon of asymmetric nuclear matter in /^-equilibrium, 
we calculate the nuclear contribution P nuc i to the total pressure of stellar mat- 
ter as well as the mass density p nuc i- Then the total pressure and total mass 
density can be easily calculated by adding the leptonic contributions. For more 
details, see ref.|31. 



3 Neutron star structure 

Neutron stars are objects of highly compressed matter so that the geometry 
of space-time is changed considerably from flat space. Thus Einstein's general 
theory of relativity must be applied. Therein the Einstein curvature tensor 
is coupled to the energy-momentum density tensor T^ v of matter (G denotes 
the gravitational constant): 

Gp, = &*GTp,{p,P{p)) (4) 

The knowledge of the EOS, i.e. pressure P as function of the energy density p 
is therefore required in order to solve eq.(|4]). For a spherically symmetric and 
static star, Einstein's equations reduce to the familiar Tolman-Oppenheimer- 
Volkoff (TOV)|U equations 



dPjr) _ Gm{r)p{r) (1 + j^yX 1 + ^mjr) ) 
dr ~ r 2 (I — 2Gm ( r ) \ ^ ' 

d ^P- = Wp(r) (6) 
dr 

For a given EOS i.e. P(p), one can solve the TOV equations by integrat- 
ing them for a given central energy density p, from the star's center to the 



star's radius defined by P(R S ) = 0. This gives the stellar radius R and the 
gravitational mass is then 

M G = m(R) =Att dr r 2 p(r). (7) 
Jo 

The case of rotating stars is more complicated, since changes occur in the 
pressure, energy density and baryon number density because of the rotation. 
In this work we adopt the method developed in ref. ||, which is a redefined 
version of Hartle's perturbative method J7| for the investigation of the gen- 
eral relativistic Kepler frequency of a rotating neutron star. We recall that 
for rotation at frequencies beyond the Kepler value Qk, mass shedding at the 
equator sets in which makes the star unstable. Therefore Qk sets an abso- 
lute upper bound on the rotational frequency. For every nuclear EOS there 
are uniquely determined values of VL K for each star in the sequence up to the 
limiting mass value. In Fig. 3 we show the results obtained with our BHF 
plus TBF equation of state for asymmetric nuclear matter and with the Dirac- 
Brueckner one. We display the gravitational mass Mq, in units of the solar 
mass M Q (M & = 1.99 10 33 g), as a function of the radius R and the central 
energy density p (in units of po = 140 MeV/ fm 3 ). The upper lying curves 
show the increase of mass due to rotation at the (absolute limiting) Kepler 
frequency, i.e. Q = VL K . The solid (dashed) line indicates a sequence of star 
models obtained with the BHF+TBF (DBHF) equation of state. We observe 
larger gravitational masses, i.e. up to 14% mass increase for BHF+TBF and 
15% for DBHF, relative to the spherical (non-rotating) Oppenheimer-Volkoff 
star model of the same p value. The large mass increase obtained is accompa- 
nied by relatively large (equatorial) radius value. Moreover one sees that, for 
a fixed value of the gravitational mass, the central star density p decreases for 
increasing values of the rotational frequency. This is because of the centrifu- 
gal force acting on the star's matter together with the nuclear force. Their 
intensity must be counterbalanced by the attractive gravitational forces. More 



details on rotating stars will be given in a forthcoming paper [16] 



4 Conclusions 

In conclusion, we computed some properties of NS's on the basis of a micro- 
scopic EOS obtained in the framework of BHF many-body theory with two 
plus three-body nuclear interactions. Our EOS with three-body forces is able 
to reproduce the correct saturation point of nuclear matter. The comparison 
with the DBHF method shows that the relativistic effects represent a partic- 
ular repulsive three-body force and gives rise to a stiffer EOS. Therefore the 



predicted values for the limiting mass and radius of neutron stars are higher in 
the DBHF case than in the non-relativistic BHF calculation. The inclusion of 
rotation substantially changes the limiting values of the gravitational mass by 
about 14-15 % for configurations rotating at their respective general relativis- 
ts Kepler frequency. A more detailed comparison between the predictions of 
non-relativistic BHF microscopic EOS and observational data on pulsars, i.e. 
fast rotation and large enough neutron star masses, is currently in progress. 
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Figure captions 



Fig.l : The energy per baryon E/A is plotted vs. the number density 
n for symmetric matter (lower curve) and for neutron matter (upper curve). 
The solid line represents a Brueckner calculation with Avl4 potential with the 
continuous choice. The squares are the solutions of the Bethe-Fadeev equations 
with the gap choice. 

Fig.2 : The energy per baryon E/A is plotted vs. the number density n in 
panel (a) for symmetric matter (lower curves) and for neutron matter (upper 
curves) . The solid line represents non-relativistic BHF calculations with three- 
body forces and the dashed line a relativistic Dirac-Brueckner one. Panel (b): 
as in panel (a) but for /3-stable nuclear matter. 

Fig.3 : The gravitational mass Mq, expressed in units of the solar mass 
Mq, is displayed vs. radius R (panel (a)) and the central energy density p 
(in units of p = UOMeV/ fm 3 ) (panel (b)) for sequences of star models 
constructed for BHF+TBF (solid line) and DBHF (dashed line) equations of 
state. The two (lower) upper lying curves refer to (non-) rotating star models. 



